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Abstract 

We present a systematical study of two-body and three-body charmless baryonic B decays. 
Branching ratios for two-body modes are in general very small, typically less than 10“®, except 
that B{B~ —> pA ) ~ 1 X 10“®. In general, B —> NA > R —> NN due to the large coupling 
constant for S;, —> BA. For three-body modes we focus on octet baryon final states. The leading 
three-dominated modes are ^ ^ pnTr~{p~), npTT^{p'^) with a branching ratio of order 3 x 10“® 
for B^ —> pn7r“ and 8 x 10“® for B^ —> pnp~. The penguin-dominated decays with strangeness in 
the meson, e.g., B~ —> ppK~^*^ and B^ —> pnK~^*\ nnK^^*\ have appreciable rates and the NN 
mass spectrum peaks at low mass. The penguin-dominated modes containing a strange baryon, 
e.g., B^ —> S'^pvT'*', S“n7r'’“, have branching ratios of order (1 ~ 4) x 10“®. In contrast, the 
decay rate of ^ —> ApTr'*' is smaller. We explain why some of charmless three-body final states 
in which baryon-antibaryon pair production is accompanied by a meson have a larger rate than 
their two-body counterparts: either the pole diagrams for the former have an anti-triplet bottom 
baryon intermediate state, which has a large coupling to the B meson and the nucleon, or they are 
dominated by the factorizable external lU-emission process. 
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I. INTRODUCTION 


Inspired by the claim of the observation of the decay modes ppir^ and ppyr+vr" in B decays 
by ARGUS in the late 1980s, baryonic B decays were studied extensively around the early 
1990s PHfBH with the focus on the tree-dominated two-body decay modes, e.g. the charmful 
decays B —>■ A^A”, and charmless ones B pp, AA. Up to now, none of the two-body 

baryonic B decays have been observed Many of the earlier model predictions are too 

large compared to experiment. For example, the previous limit on ^ pp < 7 x 10“® set 
by CLEO has been recently pushed down to the level of 1.6 x 10“® by Belle [|T^, whereas 
the model predictions are either too large or marginally comparable to the experimental 
limit (see Table II below). 

The penguin-induced charmless baryonic B decays such as B —y Up, SA have been 
studied by Chernyak and Zhitnitsky 0] based on the QCD sum rule analysis. They obtained 
the branching ratios of order (0.3 — 1.0) x 10“®. Experimentally, only the upper limits on 
B~ —>■ Ap, Ap7i~^7i~, A^p, pA (A being the antiparticle of A^^) and —>■ ApTr"*" have 
been set. 

As pointed out by Dunietz |T^ and by Hon and Soni |^, the smallness of the two-body 
baryonic decay B B 1 B 2 has to do with a straightforward Dalitz plot analysis (see Sec. 
IV for a detailed discussion) or with the large energy release. Hon and Soni conjectured 
that in order to have larger baryonic B decays, one has to reduce the energy release and 
at the same time allow for baryonic ingredients to be present in the hnal state. Under this 
argument, the three-body decay, for example B ppn, will dominate over the two-body 
mode B ^ pn since the ejected p meson in the former decay carries away much energies 
and the conhguration is more favorable for baryon production because of reduced energy 
release compared to the latter [^]. This is in contrast to the mesonic B decays where 
the two-body decay rate is generally comparable to the three-body one. The large rate of 

—>• D*~pn and B^ D*~ppTi^ observed by CLEO indicates that the decays B — 
baryons receive comparable contributions from B 
advocated by Dunietz 0. 
out recently by [^. In |21 
B- 


AppX and B DNN'X, as originally 
A theoretical study of the decay B —>■ D*pn has been carried 
we have shown explicitly that the three-body charmful decay 

■ Acp as seen experimentally . 


Acp7r~{p~) has indeed a magnitude larger than B 
By the same token, it is expected that for charmless baryonic B decays, B - 
are the dominant modes induced by tree operators and B —>• {ti, p)Bps)B 2 , B 


7r,p)B,B2 
K^*^BiB2 

are the leading modes induced by penguin diagrams. The recent hrst observation of the 
penguin-dominated charmless baryonic decay B~ ppK~ by Belle clearly indicates 
that it has a much larger rate than the two-body counterpart > pp. Of course, this does 
not necessarily imply that the three-body hnal state B 1 B 2 M always has a branching ratio 
larger than the two-body one BiB 2 - We shall examine under what circumstance that the 
above argument holds. 

In the present paper we will give a systematical study of two-body and three-body charm- 
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less baryonic B decays. The factorizable hT-exchange or hh-annihilation contribution to 
two-body decay modes is very small and hence negligible. For nonfactorizable contributions 
to two-body hnal states, we will calculate the corresponding pole diagrams at the hadron 
level. We will apply the bag model to evaluate the baryon-baryon matrix elements and 
hnd that the baryon-strange baryon weak transition is indeed dominated by penguin oper¬ 
ators. Branching ratios for two-body baryonic modes are found to be in general very small 
< 0(10“®) except for the decays with a A resonance in the hnal state. 

The study of three-body baryonic decays is more complicated. Though it in general 
receives factorizable contributions, some of them involve three-body matrix elements and 
hence are not ready to evaluate. Therefore, pole diagrams still play an essential role. The 
baryonic decay with a vector meson in the hnal state normally has a large rate which should 
be easily accessible by the existing B factories. 

The layout of the present paper is organized as follows. In Sec. II the issue of renor¬ 
malization scheme and scale dependence of Wilson coefficients is addressed. We then study 
charmless two-body baryonic B decays in Sec. Ill and compare our results with the literature 
and experiment. In Sec. IV some important three-body modes are analyzed. Sec. V gives 
discussions and conclusions. A short summary of the relevant baryon wave functions and 
the bag model evaluation of baryon-baryon matrix elements are presented in the Appendix. 


II. HAMILTONIAN 


The relevant ehective AB = 1 weak Hamiltonian for hadronic charmless B decays is 

Cl(/i)0“(/i) -F C2(/i)02(/i) VcbV^ Ci(/i)0)^(/i) C2(/i)02(/i) 


nefi{AB = l) = ^{VubV, 


a{' 


uq 


10 


VaV-J2c(p)Oi{p)\ +h.c., 


( 2 . 1 ) 


i=3 


where q = d,s, and 


or = 

03 ( 5 ) {Q^)v-A 1 )v'_a(v'+a)i 


07(9) = 2(^^)v-AEeg'W)v+A(v-A)l 
<?' 


O 2 v-A v-A 5 

O 2 {dabf3)Y_j^{q/3Ca)y_j^, 

04 ( 6 ) {Mabfi)v-A ^Z(^/3*?a)v-A(v'+A)i (2-2) 

q' 

3 

08(10) (j) V - A V4-4 (v- 4)1 


with O 3 -O 6 being the QCD penguin operators, O 7 -O 10 the electroweak penguin operators 
and {qiq2)v±A = '^lb)d 2 - The scale dependent Wilson coefficients calculated at next- 

to-leading order are renormalization scheme dependent. We use the next-to-leading Wilson 


coefficients evaluated in the naive dimensional regularization scheme [24 
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Cl = 1.082, C 2 = —0.185, C 3 = 0.014, C 4 = —0.035, C 5 = 0.009, cq = —0.041, 

Cj/a = —0.002, Cg/a = 0.054, Cg/a = —1.292, Cio/a = 0.263, Cg = —0.143, (2.3) 

at /j, = rnb(mb) = 4.40 GeV for = 225 MeV taken from Table XXII of with 
a being an electromagnetic fine-structnre conpling constant. In order to ensnre that the 
physical amplitnde is renormalization scale and 75 -scheme independent, we inclnde vertex 
and pengnin corrections to hadronic matrix elements of fonr-qnark operators p5|j26[| . This 


amonnts to modifying Ci(/i) —*• and 

T.c-(p)iQiW) = Y.‘:f{QihlA, ( 2 . 4 ) 

where the snbscript VIA means that the hadronic matrix element is evaluated under the 
vacuum insertion approximation. Numerical results for are shown in Table I (for details, 
see 1^). It should be stressed that c|® are renormalization scale and scheme independent. 
For the mesonic decay B —>■ Mi M 2 with two mesons in the hnal state, two of the four quarks 
involving in the vertex diagrams will form an ejected meson. In this case, it is necessary to 
take into account the convolution with the ejected meson wave function. 

The penguin matrix element of scalar and pseudoscalar densities, (i 3 ii 32 |^i(l ± 75 )q' 2 | 0 ), 
is usually evaluated by applying the equation of motion and it is renormalization scale and 
scheme dependent. Since the factorization scale is set at fif = rrib to obtain the effective 
Wilson coefficients listed in Table I, we will therefore evaluate the penguin matrix element 
of scalar and pseudoscalar densities at the rrib scale. 


TABLE I. Numerical values of the effective Wilson coefficients c|® for b ^ s, b ^ d and 
b ^ d transitions evaluated at pj = nib and taken from Table I of [26], where use of 

|Wb/Fcfe| = 0.085 has been made. The numerical results are insensitive to the unitarity angle 7 . 



b ^ s, b ^ s 

b^d 

b^d 

efF 

A 

1.169 

1.169 

1.169 

„eff 

-0.367 

-0.367 

-0.367 

„eff 

C3 

0.0227 -b i0.0045 

0.0226 -b i0.0038 

0.0230 -b ffi.0051 

„eff 

C4 

-0.0463 - zO.0136 

-0.0460 - i0.0114 

-0.0470 - zO.0154 

„eff 

C5 

0.0134 -b i0.0045 

0.0133 -b zO.0038 

0.0137-b ffi.0051 

„eff 

^6 

-0.0600 - i0.0136 

-0.0597 - i0.0114 

-0.0608 - i0.0154 


-0.0309 - i0.0367 

-0.0305 - i0.0324 

-0.0326 - i0.0403 

cf/a 

0.070 

0.070 

0.070 

cf/a 

-1.428 -i0.0367 

-1.428 -iO.0324 

-1.430-ffi.0403 

cfff/a 

0.48 

0.48 

0.48 


For quark mixing matrix elements, we will use |14fe/Wfe| = 0.085 and the unitary angle 
7 = 60°. In terms of the Wolfenstein parameters A = 0.815 and A = 0.2205 we have 
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rj = 0.385 cos 7 , 


p = 0.385 sin 7 , 


(2.5) 


where p and rj are the parameters in the Wolfenstein parametrization n^m 
matrix. 






III. CHARMLESS TWO-BODY BARYONIC DECAYS 


The charmless B decays into two light baryons can be classihed into two categories: the 
ones induced by the 6 —>■ m tree transition, and the ones by the 5 —>■ s penguin transition. 
The decay amplitude of i? ^ has the form 

A{B ^ B 1 B 2 ) = ui{A + B'y5)v2, (3.1) 


where A and B correspond to p-wave parity-violating (PV) and s-wave parity-conserving 
(PC) amplitudes, respectively. The decay rate is given by 


P(5^Si(l/2+)i32(l/2+)) 


+ mi + m2Ypl 

47r\ {El + mi){E2 + m2)m\ 

2 [{El + mi){E 2 + m2) +pg]^ | 
{El + mi)(^2 + m2)m| J ’ 


(3.2) 


where Pc is the c.m. momentum, Ei and mj are the energy and mass of the baryon i3j, 
respectively. For the decay B —> Bi{^~^)B 2 {^~^) with a spin-| baryon in the hnal state, the 
general amplitude reads 


A{B Bi{pi)B2{p2)) = iqf,u1{pi){C + D'y5)v2{p2), (3.3) 


where is the Rarita-Schwinger vector spinor for a spin-| particle, q = pi — P 2 and C, D 
correspond to parity-violating p-wave and parity-conserving d-wave amplitudes, respectively. 
The corresponding decay rate is 


P(R ^ Si(3/2+)S2(1/2+)) 


pI 1 i^,2 [iEi+mi){E2 + m2)+plf 

Qirmll {El + mi){E2 + m2)m% 


^,2 (me + mi + m2fpl \ 
{El + mi){E2 +1712)171% )' 


(3.4) 


As shown in Fig. 1, the quark diagrams for two-body baryonic B decays consist of 
internal VF-emission diagram, b d{s) penguin transition, hF-exchange for the neutral 
B meson and IF-annihilation for the charged B. Just as mesonic B decays, JF-exchange 
and JF-annihilation are expected to be helicity suppressed and the former is furthermore 
subject to color suppression.^] In the language of the pole model, the M^*^BiB 2 form factor 


*In contrast, IF-exchange plays an essential role in nonleptonic decays of baryons as it is no longer 
subject to color and helicity suppression. 
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is expected to be largely suppressed at = m%. As estimated by 


, the hh-exchange 


or Vh-annihilation contribution is very insignificant and hence can be neglected. The tree- 
dominated decays, e.g. B —> pp, pA~ are mainly induced by the internal hh-emission via 
b ^ u transition, while penguin-dominated modes, e.g. B~ —> Ap, TPp proceed through b —> 
s penguin transition. These amplitudes are nonfactorizable and thus very difficult to evaluate 
directly. This is the case in particular for baryons, which being made out of three quarks, in 
contrast to two quarks for mesons, bring along several essential complications. In order to 
circumvent this difficulty, it is customary to assume that the decay amplitude at the hadron 
level is dominated by the pole diagrams with low-lying one-particle intermediate states. More 
precisely, PC and PV amplitudes are dominated by ground-state intermediate states 
and I low-lying baryon resonances, respectively |^.[| This pole model has been applied 
successfully to nonleptonic decays of hyperons and charmed baryons [^,^. In general, the 
pole diagram leads to 


Bt 


mi - me* 


Bb 


9Bb^BB2 ^BbBi 
mi - ruB^ 


where 


{Bi\H^Q\Bb) — UB^aB^B^UBbi 


{Bi\H^q\BI) — iubJjbiBiUbi 


(3.5) 


(3.6) 


are PC and PV matrix elements, respectively. 

Since the weak transition does not involve momentum transfer, it can be evaluated using 
the quark model. Conventionally, baryon matrix elements are evaluated using the bag model 
or the harmonic oscillator model. In the present work we prefer to employ the MIT bag model 
H to compute the weak baryon-baryon transition for several reasons. First, it has been 
applied successfully to describe the p-wave amplitudes of hyperon nonleptonic decays and it 
is much simpler than the harmonic oscillator model for computing the PC matrix elements. 
Second, it is relatively easy to incorporate penguin operators in calculations. Third and most 
importantly, the bag model calculation gives predictions consistent with experiment, whereas 
the calculated results based on the harmonic-oscillator model are too large compared to the 
data. This will be clearly demonstrated below when we discuss B ^ pp and pA. 

However, it is known that the bag model is considerably less successful for describing 
the physical non-charm and non-bottom ^ resonances not mentioning the charm or 


bottom 2 baryon states. Therefore, we will not evaluate the PV matrix element and 
the strong coupling gB*^BB 2 S'® their calculations in the bag model are much more involved 
and are far more uncertain than the PC case [^|. Fortunately, there are some decay modes 
that are purely parity-conserving within the framework of the quark-pair-creation model 
to be mentioned shortly. Examples are B~ —> np and B NA, which will be discussed 
below. 


tphe s-channel meson pole states correspond to weak annihilation diagrams [see Fig. 1(b)]. 
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FIG. 1. Quark and pole diagrams for B ^ B 1 B 2 where the symbol • denotes the weak vertex. 
Fig. 1(a) corresponds to the nonfactorizable internal W emission or the b d{s) penguin transition, 
while Fig. 1(b) to the kF-exchange contribution for the neutral B or kF-annihilation for the charged 
B, or penguin-induced weak annihilation. 
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For strong couplings we will follow to adopt the ^Pq quark-pair-creation model in 

which the qq pair is created from the vacuum with vacuum quantum numbers ^Pq. We shall 
apply this model to estimate the relative strong coupling strength and choose = 5 

as a benchmarked value for the absolute coupling strength (see below). Presumably, the ^Pq 
model works in the nonperturbative low energy regime. In contrast, in the perturbative high 
energy region where perturbative QCD is applicable, it is expected that the quark pair is 
created perturbatively via one gluon exchange with one-gluon quantum numbers ^Si. Since 
the light baryons produced in two-body baryonic B decays are very energetic, it appears that 
the ^Si model may be more relevant. However, in the present paper we adopt the ^Pq model 
for quark pair creation for the following two reasons. First, it is much simpler to estimate the 
relative strong coupling strength in the ^Pq model rather than in the ^Si model where hard 
gluons arise from four different quark legs and generally involve infrared problems. Second, 
this model is presumably reliable for estimating the BbBB coupling when all particles are 
on their mass shell. Of course, the intermediate pole state Bb in the two-body baryonic 
decay is far from its mass shell (but not quite so in the three-body decay). In principle, 
one can treat the intermediate state as an on-shell particle and then assume that off-shell 
effects of the pole can be parametrized in terms of form factors. Such form factors are 
basically unknown, though they are expected to become smaller as the intermediate state is 
more away from its mass shell due to less overlap of initial and final hadron wave functions. 
Since we are interested in the relative strength of strong couplings rather than the absolute 
strength, it seems plausible to assume that the relative coupling strengths are essentially 
not affected by the off-shell extrapolation; that is, the strong form factors are assumed to 
be universal. We then use the experimental result for B~ —>■ Acp7^~ to fix the absolute 
coupling strength of or ||^.|3 In the future, it is important to carry out the 


more sophisticated pQCD analysis to gain a better understanding of the underlying decay 
mechanism for baryonic B decays. 

At this point, we would like to stress that although we employ the same pole-model 
framework as Jarh et al. to discuss baryonic B decays, the calculational detail is different. 
While Jarh et al. evaluated baryon matrix elements at large momentum transfer and strong 
couplings at small transfer, we consider weak transition at zero transfer and strong couplings 
at large momentum transfer as elaborated before. Another difference is related to the quark 
model evaluation of baryon matrix elements: We employ the bag model rather than the 
harmonic oscillator model. 


l-The nonresonant decay B~ Acpvr" receives main contributions from Figs. 2(a) and 2(d) shown 
in Sec. IV (see [^). In the pole model, the contribution of the former is governed by the pole. 
Therefore, a measurement of the decay rate of this mode enables us to determine the off-shell 
coupling . For detail, see f^ . 

b ^ ^ 



















TABLE II. Predictions of the branching ratios for some charmless two-body baryonic B decays 
classified into two categories: tree-dominated and pengnin-dominated. In this work, some branching 
ratios denoted by “f” are calculated only for the parity-conserving part. For comparison some 
other predictions in the literatnre are also shown. We have normalized the branching ratios to 
= 0.085 . The predictions given in [11] are carried out in two different qnark-pair-creation 
models: local and nonlocal. The line separates tree- and penguin-dominated charmless baryonic B 
decays and experimental limits are taken from [14,15]. 


Ref. § Ref. 0 Ref. |0| 


Ref. [11] 


non-local local 


This work 


Expt. 


R - 

-> pp 

4.2 X 10“® 

1.2 X 10“® 

7.0 

X 

10- 

-6 

2.9 

X 

10- 

-6 

2.7 

X 10- 

-5 

1.1 

X 

10- 

7t 

< 

1.2 

X 

10- 

-6 

B^ - 

-A nn 


3.5 X lO-"^ 

7.0 

X 

10’ 

-6 

2.9 

X 

10- 

-6 

2.7 

X 10- 

-5 

1.2 

X 

10“ 

7t 






B- 

-> np 


6.9 X 10“^ 

1.7 

X 

10’ 

-5 


0 



0 


5.0 

X 

10' 

-7 






# - 

AA 



2 X 10“ 
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ot 


< 

1.0 

X 

10- 

-6 

B- 

pA 

1.5 X lO-'^ 

2.9 X 10“'^ 

3.2 

X 

10- 

-4 

2.4 

X 

10- 

-6 

8.7 

X 10- 

-6 

1.4 

X 

10' 

-6 

< 

1.5 

X 

10- 

-4 

# - 

pA~ 


7 X 10“® 

1.0 

X 

10’ 

-4 

1.0 

X 

10- 

-6 

4.0 

X 10- 

-6 

4.3 

X 

10' 

-7 






B- 

-> nA~ 



1.1 

X 

10' 

-4 

2.7 

X 

10- 

-7 

1 X 10“ 

7 

4.6 

X 

10' 

-7 






# - 

-> nA^ 



1.0 

X 

10' 

-4 

1.0 

X 

10- 

-6 

4.0 

X 10- 

-6 

4.3 

X 

10' 

-7 






B- 

-> Ap 


< 3 X 10-6 












2.2 

X 

10“ 

7t 

< 

2.2 

X 

10- 

-6 


# 

B- 

B- 

B^ 

B- 


An 

S+p 

> SOp 

> s+A”- 
S+A- 

> S-Ao 


6 X 10“® 
3 X 10“® 
6 X 10“® 

6 X 10“® 


2.1 X 10“'^t 

1.8 X 10“®t 

5.8 X 10“®t 
2.0 X 10“'^ 


For the reader’s convenience, in Table II we give a summary of the calculational results 
presented in Sections III.A and III.B below. For comparison, some other predictions in the 
literature are shown in the same table. 

A. Tree-dominated two-body decays 

1 . pp 

As discussed before, we can neglect IF-exchange contributions to ^ > pp and simply 

focus on the internal LF-emission which is manifested as the pole diagram at the hadron level 
with the low-lying intermediate states [see Fig.1(a)]. The PV and PC wave amplitudes 
read 
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A 


9s+*^b°p 


TUr, 


ms* 


B 


rrip - msfe 


(3.7) 


Neglecting penguin contributions to the matrix element due to the smallness of penguin 
coefficients, we have 


Gf 


tP = ^yutv:Acf-cf){p\onK) 


,efr\ 


me I 


P 


x /2 


(3.8) 


for the PC matrix element, where 0i = {uh)y_^{du)y_^ and use has been made of 
(p|02|S^) = — (p|Oi|S^). The latter relation holds because the combination of the four- 
quark operators 0i -|- O 2 is symmetric in color indices (more precisely, it is a color sextet) 
and hence it does not contribute to the baryon-baryon matrix element since the baryon-color 
wave function is totally antisymmetric. In contrast, the operator Oi —O 2 is a color antitriplet 
and has isospin I = \ because the diquark ud is isoscalar due to anti-symmetrization. The 
latter feature will lead to some AJ = | rule relations, for example (|3.26|) below. 

We shall employ the MIT bag model [^] to evaluate the baryon matrix elements (see 

element 


for the method). From the Appendix of [^] or il we obtain the PC matrix 


{p|Op|E+) = -6A'(4;r), 


(3.9) 


where 


rR 


X= r^dr[uu{r)ub{r) -1- Vu{r)vb{r)][ud{r)uu{r) + Vd{r)vu{r)], 
Jo 


(3.10) 


is a four-quark overlap bag integral and Uq{r), Vq{r) are the large and small components of 


the quark wave functions in the ground (lS' 1 / 2 ) state (see, for example, |^). As stressed 
in passing, we will not evaluate the PV matrix element as its calculation in the bag 
model is much more involved and considerably less reliable than the PC one. (However, see 
2^,^] for the evaluation of PV matrix elements in charmed baryon decays.) Numerically, 


we obtain 


X = 1.52 X IQ-^GeVl 


Collecting everything together leads to 

B(B^ —>• pp)pc = 1-1 X 10“^ 

and hence 

B(B° pp) < 2.2 X 10“^ 






(3.11) 


(3.12) 


(3.13) 
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where the upper limit corresponds to Fpv/rpc = 1. Therefore, the above result is consistent 
with the experimental limit 1.6 x 10“® IIT 


We have chosen l^f 


p' 


= 5 as a benchmarked value for the strong coupling for two 


reasons. First, a calculation based on the quark-pair-creation model yields a value of 

6 ~ 10 for this coupling [^. Second, we have computed the decay B~ —>• AcP7r“ in [^. A 
£t to the measured branching ratio for this mode implies a strong coupling with the 

strength in the vicinity of order 16. Using the relation \gi^^^B-p\ = 3y^3/2 \gY,+^-B°p\ derived 
from the ^Pq quark-pair-creation model, it follows that | ~ 4.4, which is close to the 

above-mentioned model estimate. 

Note that a similar pole model calculation by Jarfi et al. yields a branching ratio of 
order 7.0 x 10“® after scaling their original result (see Table I of |T^) to \VublVcb\ = 0.085 
and to the current world average of B lifetimes [^. Since Tpv/Tpc = 0.79 is obtained by 
the same authors, and a strong coupling — 10 is used by them, it follows that 


B{W pp)pc°- = 1.0 X 10 


1-6 




(3.14) 


is predicted by Jarh et al. using the harmonic oscillator wave functions for baryons. 
Evidently, the estimate of the PC matrix element harmonic oscillator model is 

about three times as big as the one calculated in the bag model.| 


2 . K 


nn, B 


np 


The relevant intermediate states in the pole diagrams for the decays B nn and 
B~ —>■ np are A[*^ and Consider the former decay first. The PV and PC wave 

amplitudes read 


A = 
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rUr. 


rriA* 


rur, 




rUr. 




Applying the bag model leads to the PC matrix elements 

= 3y2X(47r), {n\0^^\Ab) = v^X(47r). 

For strong couplings, the ^Pq quark-pair-creation model implies |[3l[ 




(124)$^ (35) |$^T (123)$.a,(45)) 


B n 


{<I>„, (124)<*j«(35)|<*j.,(123)®v«(45)) ’ 


(3.15) 


(3.16) 


(3.17) 


§It is not clear to us how to make a direct comparison of our result for which has a dimension 

of mass, with the numerical value of a^+ shown in Table II of 11C 1 which seems to be dimensionless. 

2 ., p - 
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where the $’s are the spin-flavor wave functions and the vacuum wave function has the 
expression 


*hvac = -^{uu + dd + ss) ® + it)- 

Using the baryon wave functions given in Eq. (^I|) and the B meson wave function 


^-0 =bd®^{U - it), 


we obtain 




Consequently, 

B(B^ ^ nn) = 
Likewise, for B~ —>■ np we have 

B{B- ^ np) = gyo^-^ 
where use has been made of 

9hh^B~p . 9^0 

b 

Using the relations 




rrin - rriA. rUn - mj:. 




(3.18) 


(3.19) 


(3.20) 


(3.21) 


^\mp- tuaj ' rup - msf, 




and 






(3.22) 


(3.23) 


(3.24) 


(3.25) 


derived from Eqs. ( |3.9| ) and (|3.16|) , we And that B —>• NN amplitudes satisfy the AJ =1/2 
relation BJT0[| 


A{B^ —>• pp) — A{B —i> nn) = A{B —> np). 


(3.26) 


As mentioned before, this AJ = | relation arises because the weak operator Oi — O 2 has 
isospin J = |- 

From Eqs. ( p.l6|) , (|3.21| ) and (p.22|) , it is evident that B~ —>■ np has a larger rate than 
B^ —>■ nn. In contrast, the QCD sum rule analysis in |]^ predicts that r{B^ —> pp) > 
T{B~ —>• np) > T{B^ —>• nn). Moreover, as pointed out in ^^e decay B~ —>■ np is 
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purely parity-conserving, namely, its parity-violating amplitude vanishes provided that the qq 
pair is created from the vacuum. As pointed out by Korner |^, if the quark pair is created 
perturbatively via one gluon exchange with one-gluon quantum number {^Si model), the 
neutron in B~ np will have a positive longitudinal polarization. Therefore, a polarization 
measurement of the neutron by studying its subsequent weak decay can be used to test the 
and ^Si quark-pair-creation models. 

We are ready to compute branching ratios and obtain 

B(B^ nn)pc = 1.2 x 10“^ 

B{B- np) = 5.0 X 10“^ 








3. B° ^ AA 


Let us consider the PC amplitude oi B AA. In the pole model it receives pole 
contributions from the anti-triplet and sextet 


B = -^-^ 


rriA - m^. 


niA — rns 


(3.28) 


Using the wave functions given in Appendix A, we obtain 

(A|Or|S°) = -2X(4vr), 
(A|OfC|S;°) = -2v^X(47r), 


(3.29) 


for PC matrix elements, and 




BA 


= -V35'='o_ 


*WA 


(3.30) 


for strong couplings. Then it is clear that the PC amplitude vanishes as the mass difference 
between Sf, and is negligible. That is, this decay is purely parity violating in the 
quark-pair-creation model as noticed by Korner P| and Jarh et al. 0 some time ago. As 
noted in passing, we will not compute the PV amplitude within the framework of the bag 
model. 


4. P —> pA , nA , B —> pA , nA° 

The relevant pole diagram consists of the intermediate states for pA , pA~ modes 
(A being the antiparticle of A+’'' and likewise for other A particles) and as well as 
A[*^ for nA~ , nA° hnal states. However, it is straightforward to show that, in the ^Pq quark- 
pair-creation model, the strong coupling for A;, —NA vanishes and hence the Af, pole makes 
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no contribution. Moreover, the parity-violating part vanishes in the same quark-pair-creation 
model |]5|,p!0|. Therefore, 


D{B- ^pA—) = 


9j:+^b-a++ ®e+ 


bP 


D(W ^pA-) = 


rrip - 

9s+^b°a+ 

mp - m^. 


D{B- 


nA~) = —^^ 


mn - mj:. 


_fi _ . 9y0 \0 ^S?n 

D(^^nA^)= ” 


rrin - my. 


(3.31) 


where the PC matrix elements a^+p have been evaluated before. The relative strong 

couplings are 


= -^9, 


5's+^B-A++ - '^yT.+^-^A+ 


- \[^J^ 9y0 


^Wao 


— gyo^B-A+ — 


s+^s p • 


(3.32) 


This together with the baryon matrix elements (|3.9|) and (|3.16 ) leads to the relation 




T{B- pA—) = 3T{B- nA-) = 3T{W pA”) = 3r(5 ^ nA^), 


(3.33) 


as first pointed out by Jarfi et al. [|^. In the diquark model of [^, nA~ has a rate different 
from pA“ and nA°. Hence, experimentally it is important to test the relation (|3.33|) . 

If we apply Eq. (|3.32|) and use \gy+^-^ | = 5, we will obtain \gy+^^-^++\ = 24 and 
B{B~ —>■ pA ) = 5.8 X 10“®. Because of the strong decay A —> pvr", the resonant 
contribution from A to the branching ratio of ppTi~ would be 6 x 10“®. This already 
exceeds the recent Belle measurement B{B~ —> pp'K~) = (1.9^o;9 ± 0.3) x 10“® or the upper 
limit of B{B~ —> ppn~) < 3.7 x 10“® [^. Therefore, the coupling of the A to the B meson 


and the octet baryon is smaller than what is expected from Eq. (|3.32|) probably due to the 
different off-shellness of A. Recall that the parity-conserving transition to the A corresponds 
to a L = 2 partial wave. Therefore, the off-shell suppression on the three-point coupling of 
Sf, —RA is likely to be different from that of Aj, ^ BN. For definiteness, we will choose 




m-A++ 


= 12 and obtain 


B{B- ^pA—) = 1.4 X 10“® 


9-e+^b-a++ 


12 


(3.34) 


Thus this charmless decay B —> pA can have a large branching ratio of order 10 ® owing 
to the large coupling constant gy+^^-^++. In sharp contrast, this mode is predicted to 
be only at the level of 3 x 10“^ in the QCD sum rule analysis 0 (see also Table II). The 
branching ratios of other modes can be calculated using Eq. ( p.33|) and are shown in Table 
II. Experimentally, the decay B~ —>■ pA should be readily accessible by B factories BaBar 
and Belle. 
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B. Penguin-dominated two-body decays 


1 . B- 


Ap, B 


An 


This decay receives internal hh-emission and b ^ s penguin contributions [see Fig. 1(a)]. 
As we shall see below, it is a penguin-dominated mode. The pole diagram for B~ Ap 
consists of the intermediate states and 


A = 


9A*^B-pbA*A ^Sg*A 


ruA - niA* 


mA - mjj* 

o 


^ ^ 9a,^b-p aA,A ^ f/sg^B-pasgA 


rriA - mA, 


rriA - ms, 


(3.35) 


To evaluate the hadronic matrix elements, we notice that the combinations of the operators 
02 i+i + 02 i +2 (^ = 0, ■ • •, 4) are symmetric in color indices and hence they cannot contribute 
to the baryon-baryon matrix element. From this we can write the PC matrix element 
as 

OA.A = - cf ){A|Op|A,) - V,iV‘[(cf - cf)(A|Oy|At) 


+ (cf - cf ){A|Oy|A,) + (cf - cf ){A|OP|A,) + (cf - cS){A|Oy |A,) 


aPCi 


„efr eff 


PC I 


eff effA 


aPCi 


(3.36) 


Since the bag model implies 


{A\{sb)y_^{dd)^^JAh)pc = {A\{sb)y_^{uu)y^JAh)pc, 


(3.37) 


the baryon matrix elements of O 3 and Og can be related to Oi, while matrix element of Oj 
is related to O 5 , for example. 


(AlOg^lAfe) = (A|(s6)^_^[(mm)^_^ + {dd)^_^]\Ab)pc = -2{A\0^^\Ab) 


(3.38) 


Hence, Eq. ( |3.36| ) can be recast as 
Gf' 


OAaA = 


V2 


t4i,lC,(cf - cf ) - vw:.(-2cf + 2cf - Icf + ic")] {A|Or°|At) 


-V„V’(cf - cf + Icf - icf){A|Oy|At)|, 


Likewise, the relation 


(A|(s6)^_^(dd)^^^|S°)pc = -(A|(s6)^_^(mm)^^^|S°)pc 


(3.39) 


(3.40) 


implied by the bag model leads to 
Gf' 


“^ 2 -' ^ 


v^v:,(cf - cf) - V„V,:(-lcf + |cS)] {A|OfC|E») 


r,tv;:pf-cf){A|oy|ES), 


PC|V^0\ 


(3.41) 
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Therefore, the PC matrix element for — A weak transition does not receive QCD penguin 
contributions. 

Applying Eqs. (P^) and (p5|) we obtain 


(A|Or"|A,) = -Xi(47r), 

(A|Oy|At) = ^271 + 2K2 - Y; + Yi)(i 7 ,), 


(3.42) 


(A|Or^|S^) = -^(Xi + 3X2)(4vr), 

(A|Oy|E») = -^IVi + n - 2(K,' - Yi)][Air), 


in the bag model, where 

rR 


Xi = / r^dr[us{r)vu{r) - n^(r)ti„(r)][M„(r)vb(r) - n„(r)Mf,(r)], 


/o 


rR 


X2 = / r^dr[us{r)uu{r)+Vs{r)vu{r)][uu{r)ub{r)+Vu{r)vb{r)], 


rR 


Yi= r^dr[us{r)vb{r) - Vs{r)ub{r)][uu{r)Vu{r) - n„(r)M„(r)], 


rR 


¥( = / r‘^dr[us{r)vb{r)+Vs{r)ub{r)][uu{r)vu{r)+Vu{r)uu{r)], 


rR 


¥2= ¥dr[us{r)ub{r)+Vs{r)vb{r)][uu{r)uu{r)+ Vu{r)vu{r)], 


rR 


Y 2 = / r^dr[us{r)ub{r) - Vs{r)vb{r)][uu{r)uu{r) - 


(3.43) 


are four-quark overlap bag integrals. Finally we arrive at 
B(B- ^ Ap) = -4=9 v+^tP, (^ 




rprt)y+ I ' i ’ 

V2 P \mA- niA^ tua - mT.J 


(3.44) 


where use has been made of Eqs. (|3.23|) and (|3.24|) . 

The bag integrals have the values 

= -4.6 X 10“® GeV^ X 2 = 1.7 x 10“^ GeV^ Fi = 0, 

¥; = 4.5 X 10“^ GeV^ F 2 = 1-7 x 10“^ GeV^ ¥' = 1.2 x 10“^ GeVl (3.45) 

It is easy to check that oa^a and hence the decay is penguin dominated. For the branching 
ratio we hnd 


}3{B- Ap)pc = 2.2 X 10 


-7 




(3.46) 


For B —>■ An, it has the same rate as B Ap. 
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S+p 


2 . K 


We consider the pole diagram with the intermediate states 


A = 




ms - ms* 

o 


B = 


rriY, — mj.+ 


(3.47) 


The PC weak matrix element for — E"*" transition reads 


Gf 


“e+e* - ^ 


niK:.(cf - cf) - vav:,(-cf + cf - cf + c")l (e+iopie+) 


■VaV;,(cl 


cf + c% 


cf){E+|OP|Ep . 


(3.48) 


In the bag model, 


(E+|Op|Ep = 3(A',-9Ay(4i), 
{e+| 05 '’<=|e+) = |[n + + 4(f; - n)l(4^). 


We obtain numerically 

E+p)pc = 1.8 X lO" 




(3.49) 


(3.50) 


Note that the branching ratio is predicted to be 5 x 10 ® in the QCD sum rule analysis of 
1^, which is larger than our result by two orders of magnitude (see Table II). 


3. B- 


SOp 


The intermediate low-lying pole states for this decay are and 


A = 


5'A*^B-p ^A*S0 ^S0*S0 


ms - mA* 


ms - ms* 


^ _ gA,^B-p QA^SO ^ 


ms - mA, 


ms - ms, 


The PC matrix elements are given by 
Gp ^ 


ceje. - ^ 


VMcf - cf) - Vt,v:,{-2cf + 2cf - icf + iclS)] {E“|OnE») 


-V,tV-{cf - cf + Icf - icf )(E»|OP|E»)|, 


(3.52) 


and 


_ Gf 

CA»E. - 


VutV:,(cf - cf) - V,cV;,(-^cf + |cS)] (E'lOriAt) 


r,i,v;:(cf - cf XE^ioriAi) 


OiryPCi 


(3.53) 
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where in the bag model 


(E"|Or'-|Ai> = -^(.Y, + 3A'2)(4;r), 

(E»|Or|Aj> = -^[F, + Y 2 - 2(Y; - F')|(4;r), 
(E»|OnE»> = l(A,-9A2)(47r), 

(E«|Or|E»> = |[F + F + 4(F - F')](4>r). 


(3.54) 


Hence, 


B = — 


3^3 oa, 


so 


+ 


®eOso 


P \mj: - rriA,^ uit. - J ’ 

where use of Eqs. (p.23|) and (|3.24D for strong couplings has been made. We obtain 


^T.+^Wp 


B{B- S°p)pc = 5.8 X 10" 

Again, the QCD sum rule prediction for this mode is much higher [0. 

4. B- S+A—, E-A°, S+A- 


(3.55) 


(3.56) 


As stated before, the decays B —> S+A , S A°, B 
conserving contributions |]5|Jl0| so that 

V^-I- A N fi's+^B-A++®S+S+ 

S+A”) = ^—, 


S'^A only receive parity- 


D{B- 

D(^ 

D{B- 


mj: - 

9s+^-^a+ ®s+s+ 
rriY. - ms(, 

^ ^s-^g-AO Qs-s- 


S+A-) = 


(3.57) 


rriY, - ms, 

The PC matrix element as+s+ been evaluated before and as“s- ~ ®s+s+- strong 
couplings we get 


5's+^a++b- 


^s+^A+m 


in the ^Pq model. Collecting all the results gives 
B{B- S+A”) = 2.0 X 10“^ 

B(^ S+A-) = 6.3 X 10“® 

B{B- S-A°) = 8.7 X 10“® 



- = 2/ 

^S+- 

■*A++B- 


12 

^S+- 

aA++B- 


12 

l/s+- 

^A++B- 




(3.58) 


12 


(3.59) 


where we have followed the discussion of P —> A^A to choose the coupling | = 12 

b 

as a benchmarked value. 
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C. Comparison with other models 


As discussed in passing, though we adopt the same pole-model framework as Jarh et 
al. |]ro| for describing two-body baryonic B decays, a crucial difference is that weak baryon 
matrix elements are evaluated by Jarh et al. at large momentum transfer and strong cou¬ 
plings at small transfer, whereas the weak transition is computed at zero transfer and strong 
couplings at large momentum transfer in our case. In general, the difference in numerical 
results shown in Table II comes mainly from the fact that we use the bag model rather than 
the harmonic oscillator model to evaluate weak baryon transitions. 

In the following we compare our results with the diquark model by Ball and Dosch 
||Tl| and the QCD sum rule analysis by Chernyak and Zhitnitsky (see also Table II). For 
B NN decays, the diquark model has one unique prediction, namely, there is no B~ —> np 
decay, while pp and nn hnal states have the same rates. In contrast, the sum rule approach 
predicts that r(B^ —> pp) > r{B~ —>• np) > r{B^ —>• nn) (see Table II), while in our case 
r{B~ — np) > r{B —> pp) ^r{B —>■ nn). Therefore, a measurement of the relative rates 
oi B NN (especially B~ —>■ np) will serve to test the three models. 

As for the tree-dominated modes B ^ NA, they are suppressed in the diquark model 
because the operators Oi and O 2 can only generate scalar diquarks whereas the decuplet 
baryons are made of a vector diquark and a quark. Likewise, they are also suppressed in the 
sum rule analysis. In sharp contrast, these modes have sizable branching ratios in the pole 
model, namely, r{B —>■ NA) > r{B —>■ NN), owing to the large coupling of the intermediate 
state Sb with the B meson and the A resonance. 

The penguin-dominated decays have smaller rates than B —y pp in the diquark model as 
the penguin operators are not included in the original calculations by Ball and Dosch (the 
effect of the penguin operators in this model was recently discussed in [Q). In contrast, the 
sum rule approach predicts branching ratios of order (2 — 6) x 10“® for B —> Ap, T,p, SA. 
In our work, the decay rates of penguin-dominated decays are in general small. 

In short, measurements of the relative rates of i? NN, Sp, SA will suffice to differ¬ 
entiate between above-mentioned three models. 


IV. CHARMLESS THREE-BODY BARYONIC DECAYS 


As noted in the Introduction, the study and search of the three-body baryonic B decay 
B —y B 1 B 2 M with M being a meson are mainly motivated by the experimental observation 
that B{B~ — >• Acpvr”) > B{B^ — >• Acp) [Q and B{B~ ppK~) > B(B^ pp) [p3| . 
Theoretically, it has been argued that the emitted meson M in the three-body hnal state 
carries away much energies and the conhguration is more favorable for baryon production 
because of reduced energy release compared to the latter |T^ . Roughly speaking, the reason 
that the two-body baryonic decay B B1B2 is smaller than the mesonic counterpart B — 
Ml M 2 stems from the fact that one needs an additional quark pair production in the internal 
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ly-emission diagram [Fig. 1(a)] and two qq pairs in weak annihilation diagrams [Fig. 1(b)] 
in order to form a baryon-antibaryon pair. A qq prodnction is snppressed by either a strong 
conpling when it is prodnced pertnrbatively via one glnon exchange or by intrinsic softness 
of nonpertnrbative pair creation [^. In the three-body baryonic decay, the emission of the 
meson M will carry away energies in snch a way that the invariant mass of B 1 B 2 becomes 
smaller and hence it is relatively easier to fragment into the baryon-antibaryon pair. 

One can also nnderstand the above featnre more concretely by stndying the Dalitz plot. 
Dne to the V — A natnre of the b —>■ udu process, the invariant mass of the diqnark ud peaks 
at the highest possible valnes in a Dalitz plot for b udd transition (see [^ and footnote 
[91] in [^). If the ud forms a nncleon, then the very massive udq objects will intend to 
form a highly excited baryon state snch as A and N* and will be seen as Nn'K{n > 1) ||16|| . 
This explains the non-observation of the NN hnal states and why the three-body mode 
NNti{p) is favored. Of conrse, this does not necessarily imply that the three-body hnal 
state B 1 B 2 M always has a larger rate than the two-body one B 1 B 2 . In this section we will 
stndy some leading charmless three-body baryonic B decays and see nnder what condition 
that the above argnment holds. 

The qnark diagrams and the corresponding pole diagrams for decays of B mesons to 
the baryonic hnal state B 1 B 2 M are more complicated. In general there are two external 
VF-diagrams Figs. 2(a)-2(b), fonr internal kF-emissions Figs. 2(c)-2(f), and one IF-exchange 
Fig. 2(g) for the nentral B meson and one IF-annihilation Fig. 2(h) for the charged B. 
Becanse of space limit, pengnin diagrams are not drawn in Fig. 2; they can be obtained 
from Figs. 2(c)-2(g) by replacing the b ^ u tree transition by the b —> s{d) pengnin 
transition. Under the factorization hypothesis, the relevant factorizable amplitndes are 


Figs. 2(a), 2(c) : A (x {M\{q 3 q 2 )\ 0 ){BiB 2 \{qib)\B), 

Figs. 2(b), 2(d) : Flex (SiB 2 |(gig 2 )| 0 )(M|(g 3 &)| 5 ), (4.1) 

Figs. 2(g), 2(h) ; Flex (SiB 2 M|(gig 2 )| 0 )( 0 |(g 3 &)|i?). 


Since the three-body matrix elements are basically nnknown, only the factorizable amplitndes 
for Fig. 2(b) or 2(d) are calcnlable in practice. 

The tree-dominated three-body modes of interest are: 

^ 7r+(p+){np, AS”, D'A, “”5°, pA”- 

B° ^7i-{p-){pn, S+A, AS+, A++p, ■ ■ ■}, (4.2) 

B~ —>■ vr“(p“){pp, nn, SS, AA, AA, • • 

while some interesting pengnin-dominated decays are 

^ —>• 7r^(p’'"){Ap, S°p, S“h, S“A, S^A 
5“ ^ 7r-(p-){S+p, An, S°h, S°A, S+A",---}, 

S+A, AS+, A++p, ■ • ■}, (4.3) 
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FIG. 2. Quark and pole diagrams for three-body baryonic B decay B —> B 1 B 2 M, where the sym¬ 
bol • denotes the weak vertex. Figs. 2(a) and 2(b) correspond to factorizable external IF-emission 
contributions, Figs. 2(c) and 2(d) to factorizable internal VF-emission, Figs. 2(e) and 2(f) to 
nonfactorizable internal VF-emission, Fig. 2(g) to VF-exchange and Fig. 2(h) to VF-annihilation. 
Penguin contributions are obtained from Figs. 2(c)-2(g) by replacing the 6 —> u tree transition by 
the b —> s{d) penguin transition. 
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B- ^ K-^*\pp, nn, EE, AA, AA, • • 
^Tf^*\pp, nn, EE, AA, AA, • • ■}, 

B- ^Tf^*\np, AS", E°S-, E'A, pA”,---}. 


In the present paper we will focus on octet baryon final states. 

To evaluate the factorizable amplitude for Fig. 2(b) or 2(d) we need to know the octet 
baryon form factors defined by 






mi + m 2 


(^fMuq + 


mi + m 2 


% 


± 




mi + m 2 




mi + m 2 




75[^^2(P2), (4.4) 


where q = pi+p 2 - For octet baryons one can apply SU(3) symmetry to relate the vector form 
factors to the nucleon magnetic and electric form factors. In general, SU(3) symmetry 
implies 




^fie=(g2) = (4.5) 


where and well-known Clebsch-Gordon coefficients and FY[q^) and 

DY {qY are reduced form factors. The nucleon matrix element of the electromagnetic current 
is given by 


(iV(pi)A(p2)|Jr|0) =h^(pi) FYYh, + ^ 


F2{qY 


2 m 


N 


-cTfpuq 


Vn{P2 


(4.6) 


= ^U(3) symmetry allows us to determine FY[qY and DY{qY sepa¬ 

rately. The results are (see e.g. p6| ) 


Since J,®™ 

H' 


A2(«) = 

DLm = -Y"M 

F^(t) = D'pt) = 0, 

(4.7) 

with t = q^. At f = 0 we have 




F,'"(0) = 1, 71^(0) = 0, 

F 2 (0) = + 2^^., 

G2''(0) = -A„. 

(4.8) 


where Kp = 1.79 and Kn = —1.91 are the anomalous magnetic moments of the proton and 
neutron, respectively. 

The experimental data are customarily described in terms of the electric and magnetic 
Sachs form factors GY(t) and GYiit) which are related to Ff and F^ via 


Gg"(i) = G5f (t) = FrW + (4.9) 

4m^ 

A recent phenomenological fit to the experimental data of nucleon form factors has been 
carried out in using the following parametrization: 
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irP (fW _ f , ^2 2:3 a ;4 xs 




-7 


|GmWI = (|+|) 


In^T 


n -7 


Qf 


(4.10) 


where Qo = Aqcd 
fit values 


300 MeV and 7 = 2 + ^ = 2.148. We will follow to use the best 


xi = 429.88 GeV^ xa = -10783.69 GeV®, X 3 = 109738.41 GeV®, 


X 4 = -448583.96 GeV^°, X 5 = 635695.29 GeW^ 


12 


and 


yi = 236.69 GeV^ 


j /2 = -579.51 GeV®, 


(4.11) 


(4.12) 


extracted from neutron data under the assumption |G^| = \G^\. Note that the form factors 
given by Eq. ( |4.10| ) do satisfy the constraint from perturbative QGD in the limit of large 
t [^. Also as stressed in ||^, time-like magnetic form factors are expected to behave like 
space-like magnetic form factors, i.e. real and positive for the proton, but negative for the 
neutron. 

A new empirical £t to the reanalyzed data for G\j{t) in the region 0 < t < 30 GeV^ is 
recently given in [p7|] : 

Up 


GUQl = 


1 + ziQ Z2Q‘^ -\- z^Q^ Zi^Q^ -\- z^Q^ 


(4.13) 


with 


zi = (0.116 ± 0.040) GeV\ Za = (2.874 ± 0.0.098) GeV^ Z 3 = (0.241 ± 0.107) GeV®, 
Z 4 = (1.006 ± 0.069) GeV-^ Z 5 = (0.345 ± 0.017) GeV® , (4.14) 

and Up = 2.79. An empirical fit to the proton electromagnetic form factor ratio is also 

G^t) 


presented in |37 




GW) 


= 1.0 - (0.130 ± 0.005)[« - (0.04 ± 0.09)], 


(4.15) 


for the range 0.04 < t < 5.6 GeV^, indicating that the form factor ratio decreases with 
increasing Q^. 

As for the axial form factors, no useful information can be extracted from SU(3) sym¬ 
metry. Nevertheless, perturbative QGD indicates that, in the range of high Q^, the form 
factors flit) and 5 'i(t) dominate a.t t ^ oo and all others are suppressed by powers of m/Q 
3^. Moreover, all octet-octet and octet-decuplet form factors at large t can be related to 


the magnetic form factors of the nucleon G\j and G^j (see Tables II-IV of [^). Hence, the 
axial form factor gi at large momentum transfer is fixed. 
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A. Tree-dominated three-body decays 


1. B —np7r+(p+) 

This decay receives factorizable contributions from Figs. 2(b), 2(d) with b ^ d penguin 
transition, 2 (g) and a nonfactorizable contribution from Fig. 2(e). As the two-body baryonic 
decay, we can neglect the IF-exchange contribution. Moreover, we shall assume that this 
mode is dominated by the factorizable term from Fig. 2(b) as it is governed by the parameter 

di 

A(B^ np7r+(p+))fact = ^VudV*t,ai{7r^{p'^)\{ub)^_jB^){np\{du)^_jO), (4.16) 

where ai = c'f + c|®/3 and we have neglected penguin contributions because the penguin 
Wilson coefficients C 3 , • • •, Cio are numerically very small. The two-body meson matrix ele¬ 
ments are parametrized in terms of the form factors Fq and Fi for B — tt transition 

{ir+(p,)|(t<t)„_JB"(pB)) = Ff(q^)(pB+P,)„ + (<'(«") - f,*'(«")) 
and form factors V, Aq, Ai, A 2 for B — p transition 


(p^(Pp)IWv-aI^(Pb)) =- 7 - -il {ttib + mp)e*^Af%q^) (4.18) 

TTIb “T 'dip 


^ {PB+Pp)tqA^^{q‘^) - 2mp ^ ^^ qp [A^^(g^) - A^^iq^ 

TTIq ~\~ Trip Q 


with q = Pb — Ptt{p) and 


2 x ^B+^p Bpf 2\ niB — m 


A^T) = ^^^^Af^(g^) - 




(4.19) 


The factorizable amplitude for the pion emission reads 

G 

A(^ np7r+)fact = -^VudV*h aiUn [(a;i^ + b) - + d)'^^] Vp, (4.20) 


where 


a = 2 /r(t)Ff-(t) + Af^^it)Ff-it), 

b = -2/r(t)Ff-(t)(p,, - pp) ■ pJi2m^) + - ml)/i2m^), 

c = 2 g 7 it)Ff-it), 

d = 2mr,g7{t) Ff^t) + {F^^t) - Ff-(t)) ^^~^' 

-2gr{t)F,^7t){Pn - Pp) ■ Pn/{2m^) + g7{t)Ft mml - ml)/{2m^), (4.21) 
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and t = = {j)B — V-kY = {Pn + PpY- The amplitude for the p meson case is more 

cumbersome. 

Since the relevant Clebsch-Gordon coefficients are = 1, it follows from Eqs. 

([4.5|) and (|4.7|) that the weak form factors have the form 


fS(t) = Fp_{t) + = FIM - Fl,_(t). 


(4.22) 


In terms of the nucleon magnetic and electric form factors, the weak form factors read 


/^(f) = 




4m1 




t/{Am%) - 1 


t/(4m^) - 1 


^ G%{t)-GUt) ^ G^t) - GUt) 


(4.23) 


f/(4m^)-l t/(4m^)-l 

According to perturbative QCD, the weak form factors in the large t limit have the expres¬ 
sions 1^ 


/r(() ^ g7(t) - + G’i,((). 


(4.24) 


It is easily seen that this is consistent with the large t behavior of given by Eq. ( 4.23 ). 

The total decay rate for the process B {ps) —> n(pi) +p{p 2 ) + is computed by 


r = 


1 


(27r)3 32m| 


\A\'^ dm\2dm\ 




(4.25) 


B 


where rnfj = {pi + pjY with ps = To make a numerical estimate, we apply two different 
empirical hts of G\j{t)\ Eq. (|4.10|) denoted by CHT (Chua-Hou-Tsai) and Eq. ( [4.13| ) denoted 


by BKLH (Brash-Kozlov-Li-Huber). For the proton electric form factor, we shall follow 


to assume |G^(t)| = \G\^{t)\ for CHT form factors and Eq. (|4.15|) for BKLH form factors. 
That is, we assume that Eq. (|4.15|) is applicable also to the large t region. As for B — 7r(p) 
form factors, we consider two distinct models: the Bauer-Stech-Wirbel (BSW) model 
and the Melikhov-Stech (MS) model based on the constituent quark picture |^.|^ The BSW 
model assumes a monopole behavior for all the form factors. However, this is not consistent 
with heavy quark symmetry for heavy-to-light transitions. For example, the form factors 
Fi, V, Aq, A 2 in the infinite quark mass limit should have the same dependence and they 
differ from Fq and Ai by an additional pole factor [^. Nevertheless, we apply this model 
for comparison. 


**The QCD sum rule method based on the light-cone sum rule analysis [40| is also one of the 
popular form-factor models. However, we found that some divergence occurs in the phase space 
integration when applying this model. 
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Considering only the vector-cnrrent contribntion to baryon matrix element, we obtain 


the resnlts shown in the hrst entry of Table III. Onr calculations are in agreement with 118 


when the BSW model and CHT form factors are used. However, we see from Table III that 
the branching ratio for ^ npp^ in the BSW model is slightly larger. This is ascribed to 
the monopole form factor dependence for all the B — p form factors. If one changes the 
form factor momentum dependence from monopole to dipole form for Ai and V (sometimes 
referred to as the BSWII model in the literature), the resulting branching ratios are very 
similar to that in the MS model. 


TABLE III. Branching ratios of B —> npTT^{p'^) in two different form-factor models for B—7t{p) 
transition. Two distinct empirical fits for the proton magnetic form factor given in Eqs. ( |4.10| ) 
and ( ]4.13| ), denoted by CHT and BKLH respectively, are utilized. The neutron form factors are 
taken from ( 4.10| ) and ( 4.12| ). Branching ratios in the first entry are without contributions from 
the axial form factors and those in the second entry take into account contributions from 


the asymptotic form factor given by Eq. (4.24). 



'^M 

(CHT) 

G^m 

(BKLH) 

MS 

BSW 

MS 

BSW 

B - 

-> npir^ 

1.7 X 10“® 

1.8 X 10“® 

8.0 X 10“^ 

8.5 X 10“^ 



1.7 X 10“® 

1.9 X 10“® 

1.1 X 10“® 

1.3 X 10“® 

b” - 

-> npp^ 

3.3 X 10“® 

4.8 X 10“® 

4.2 X 10“® 

5.5 X 10“® 



3.4 X 10“® 

4.9 X 10“® 

4.6 X 10“® 

5.9 X 10“® 


To estimate the contribution from the axial vector current, we might assume that gi{t) 
takes the asymptotic form |G^(t)-|-G')(^(f) [see Eq. ( |4.24|) ]. It turns out that the contribution 
due to gi{t) is very small for the CHT form factor but not negligible for G^(BKLH). It 
is interesting to notice that the rate of npp'^ is larger than that of npn'^ by a factor of 2 ~ 3 
if the CHT parametrization for is employed, whereas the ratio becomes as large as 5 for 
G^(BKLH). 

Since both and B can decay into np7T~^{p'^), experimentally one has to disentangle 
the “background” contribution from the B^ — B^ mixing or to tag the B meson. Therefore, 
we will give an estimate of —>■ pfni~{p~) next. 


2. H —>• pnvr (p ) 


This decay receives contributions from Figs. 2(a), 2(e) and 2(g). As the previous de¬ 
cay, we will assume that it is dominated by the factorizable contribution from Fig. 2(a). 


Unfortunately, as shown in Eq. (^), it involves a three-body matrix element that cannot 
be evaluated directly. Instead, we will evaluate the low-lying pole diagrams with the strong 
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process B 
1 
2 


followed by the weak decays {Aj,*-*, ^ tt Consider 

the intermediate poles. Applying factorization to A^ 


(*) 


A(^ 


pmi 


— - T/* 

~ ^^^^udYub9 


A "r^ 

Ai)^B n 


vr p yields the pole amplitnde 

fn cii Up\ hp^ ■ Pp + - mp 


75 


+g^‘’^{ml) 2p^ ■ Pp - Ai'rnAh + rrip) 


>Vn X 


{Pp + P^y 


(4.26) 


where we have employed the heavy-light baryon form factors dehned by 
{piPp)\{ub)y^JAb{pAj) = Up\ /^^(p^)7/. + 


+ rrip 


rriA^ + rup 


Pirp 


± 


A,pi 2 ^ , • 

91 [pjip + — rz:r^^^>'P^ +-^- P^t^ 


rriA. + mp 


rriA, + rrip 


75^«a„ (4.27) 


with p^ = Pa, - Pp. 

For the heavy-light form factors and we will follow [^] to apply the non- 

relativistic qnark model to evalnate the weak cnrrent-induced baryon form factors at zero 
recoil in the rest frame of the heavy parent baryon, where the quark model is most trustwor¬ 
thy. This quark model approach has the merit that it is applicable to heavy-to-heavy and 


heavy-to-light baryonic transitions at maximum Following we havef^ 


= 97(9™) = 0.86, ffiilPl = 93“’(9™) = -0.51, 

= aHii) = -0.22, 

for Ah — p transition at zero recoil = (uta^ — mpY, and 

ffyal) = 1-65, fpyi,) = 1.92, f?yi) = -1.72, 
gF)?™) =-0.17, 9?’’(9™) = 0.04, 9““’’(9."J=0.10, 


(4.28) 


(4.29) 


for — p transition at g^ = {mj:, — m^)^. Since the calculation for the g^ dependence 
of form factors is beyond the scope of the non-relativistic quark model, we will follow the 
conventional practice to assume a pole dominance for the form-factor g^ behavior: 


bThere is another pole diagram with the weak decay B — > 7r“7r'’'(p'’') followed by the strong 
process 7r^(p''“) —> pfi [see Fig. 2(a)]. However, this pole amplitude is expected to be suppressed 
as the intermediate pion state is far off its mass shell. Consequently, the nucleon-nucleon-pion 
coupling is subject to a large suppression due to the form-factor effects at large . 


b-The form factors for Ah — p transition at g^ = 0 are given in Table I of |44|. For — p form 
factors at zero recoil, it can be evaluated using Eq. (22) of |Q. Note that the spin factor is p = — | 
and the flavor factor is A^s^’p = l/\/6 for — p transition. 
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fiq^) = f{ql) 


1 — /m\ 


(4.30) 


where my (m/i) is the pole mass of the vector (axial-vector) meson with the same quantum 
number as the current under consideration. The function 


G{q^) = 


'1 


pole 


1 -9V™pol. 


(4.31) 


plays the role of the baryon Isgur-Wise function C(<^) ^Q' transition, namely, 

G = 1 at = q^. Previous model calculations of C(<^) P5|-|i5[] indicates that it is consistent 
with G{q^) with n = 2. However, a recent calculation of C('^) 01 yields 

1.23+0.4/tj 


m = ( 


1 -|- o; 


(4.32) 

ppK~ 


and this clearly favors n = 1. As we shall below, the recent hrst observation of B~ 
by Belle also favors a monopole dependence for baryon form factors. 

The calculation of ^ > pnp~ is similar to that of pn7i~ except that the vacuum-p matrix 

element now reads 


{p \d-ff,u\0) = fprUpE*^, (4.33) 

and that the computation is much more tedious than the pion case, though it is straight¬ 
forward. Using the pole masses my = 5.32 GeV, m^ = 5.71 GeV and the decay constant 
fp = 216 MeV, we obtain 

13(B^ ^ pfiTT-) = 2.8 X 10-^ (2.7x10-^), 

13(B^ pnp-) = 8.2 X 10“® (8.2 x 10“^), (4.34) 


for a monopole (dipole) q^ dependence for baryon form factors. Since g 


^b- 


.*'^n 


3^/3 [cf. Eq. (p.20|) ], the contribution due to the A;, and poles is destruc¬ 


tive. In the calculation we have used \q. ^ I = 16 [21 . 

Three remarks are in order. First, in the calculation we have neglected other nonfac- 
torizable contributions from Fig. 2(e). For the pole diagrams, we did not evaluate the | 
pole contributions owing to the technical difficulties for the bag model in dealing with the 
negative-parity baryon states. Second, since n = 1 is favored by the recent measurement 
of the decay B~ —> ppK~, as we shall see below, it turns out that B 


0 


^0 


branching ratio of order 1 x 10“^ for n = 1. Third, the decay B 
nant contribution B^ —>■ pA~ followed by the strong decay A 
ratio for ^ pA~ is only of order 6 x 10” 
to the A is thus negligible. 


npp'^ has a large 
—>■ pnn~ receives the reso- 
nTT~. Since the branching 
see Table II), the resonant contribution due 


3. ^ {E°S-, E-A, 5-S°}7r+(p+) 
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The calculation for the decays B S“A, is the same as that for 

B —i> np7r+(p+) except for different baryonic form factors in the hnal states. The relevant 
Clebsch-Gordon coefficients for weak form factors are (see e.g. PS[1) 




12 

3’ 


/-'” = 0, 

Then (^ ) and lead to 


= 0 , 

= V2, 


= 1 , 


= - 1 . 


(4.35) 


/tf4«) = -\jl fTM ft?* = -K-M - 2F;:,(t), 


/& (‘) = 


(4.36) 


A straightforward calculation gives B{B^ S“A7r+) = 2.9 x 10“^, B(B^ S“S°7r+) = 
2.0 X 10“^ and B(B^ —> S°S“7r+) = 6.4 x 10“®. Compared to the npn^ mode, the decay 
rates of above three decays are suppressed owing to smaller baryon form factors and less 
three-body phase spaces available. 


4. B —>• {pp, nn, S+E (p ) 


Let us first consider the decay B —>■ ppvr . It receives factorizable contributions from 
Figs. 2(a) and 2(d): 

A{B- ^pp7r")fact = ^VubV*^[ai{Tt-\{du)y_^\Q){pp\{uh)y^^\B-) 

+ a2(7r“|(d6)^_^|5“)(pp|(hM)^_^|0)} = Ai A 2 , (4.37) 

where 01^2 = c®®-l-c|®/3. In analog to the previous mode, we will evaluate the corresponding 
low-lying pole diagrams for the factorizable external hF-emission amplitude, namely, the 
strong process B~ —{A[,*\ S°*'*^}p, followed by the weak decay —>• pvr". Its 

amplitude governed by the poles is given by 


G 

Til = -^VubV*agx^^B-pf7r ai Up{f^^^{ml)[2p^ ■ Pp + ^rriA^ - mp)]'^^ 

+ 9i'’^{ml)[2p^ ■ Pp - ^n{mA^ + mp)]}vp x - -;-^ 

yPp -h PttJ — 

+ (Afc ^ SO), (4.38) 


where we have applied factorization to the weak decay {A;,, S^} —>■ p7r . To evaluate the 
factorizable amplitude A 2 , we apply the isospin symmetry relations^ 


§§This isospin relation amounts to assuming (A|(ss)^,_^|A) = 0, an assumption supported by the 
OZI rule. 
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(4.39) 


{n\{uu)^_^\n) = {p\{dd)^_^\p), {n\{dd)^_^\n) = (p|(mm)^_^|p), 

to relate the form factors /f^ and /f^ appearing in the vector current pp matrix element 


(p(Pi)p(P2)|m7m«|0) = Up{pi) /r(g^)7M + * 




2m„ 


-CTpuq 


Vp{p2), 


(4.40) 


to the electromagnetic form factors Fi and F 2 dehned in the nucleon matrix element Eq. 
(|;D. We End 


/f>(i) = 2Ff(«) + F^{t), /f («) = 2Ft{t) + F^(t). 


(4.41) 


A straightforward calculation indicates that the contribution from 02 is small and negli¬ 
gible due mainly to the small vector form factors fi^ 2 - The ai contribution gives a branching 
ratio of order 3.8 x 10 “® for n = 1 and 2.7 x 10 “^ for n = 2 . As we shall see below, as far 
as the factorizable Oi contribution is concerned, the tree-dominated B~ —>■ pp7r~ and the 
penguin-dominated decay B~ —>■ ppK~ have almost the same rate and the latter has been 
observed recently 


In some sense this is very similar to the mesonic decays B ^ Kn 
and TTTT. Without the chiral enhancement for penguin contributions, one will have tttt > Kn. 
The experimental observation that K~ti^ > and K 7t~ > clearly implies the 

importance of penguin chiral enhancement. It is quite possible that for baryonic B decay 
we also have ppK~ > ppn~. Note that the 02 contribution to ppn~ is destructive and it 
is subject to many uncertainties. For example, the axial-vector current contribution to the 
pp matrix element has been neglected so far and the value of 02 is numerically very small if 


02 = C 2 -t-c®"/3. A large value of 02 of order 0.40 — 0.55 as indicated by the recent obser¬ 
vation of ^ I^, and an inclusion of axial form factor contributions may suppress 

pp'K~ relative to ppK~. Another effect we have neglected thus far is the penguin contribu¬ 
tion. Just as the B —>■ tttt decay, the tree-penguin interference for B~ pp7T~ may turn out 
to be destructive for a certain range of the phase angle 7 . In view of the aforementioned 
considerations, we will prefer to carry out a full analysis of B~ {pp, nn, • • •} 7 r“(p“) 

decays in a separate publication. It appears to us that B~ pp7i~ should have a branching 
ratio at least of order 10“®, based on the recent measurement of B~ —> ppK~ to be discussed 
below. 

Thus far we have focused on the nonresonant decay of ppn~. It also receives resonant 
contributions, for example B~ —>■ pA and B~ —>■ pA^‘^(1440). As discussed in Sec. III.A, 
the branching ratio of the former is of order (1 — 2 ) x 10 “®, to be compared with the recent 


measurement by Belle [23 


B{B- ppTT-) = (1.9lJ;^ ± 0.3) x 10 


1-6 


(4.42) 


Therefore, the direct nonresonant contribution is probably smaller than the resonant ones. 
Experimentally, it is thus important to study the resonance effects through the Dalitz plot 
analysis. 
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B. Penguin-dominated three-body decays 


1. B ^ NNK^*'> 


The decay B —>■ ppK is mainly governed by the diagrams Figs. 2(a) and 2(c) with 

the factorizable amplitude 


A{B —> ppK 




with 


Gf 

72 




a2{pp\{uu)^_jO){K W|(s6)^_^|5 ) 


- VttV,: 


a 3 {pp\{uu + dd +ss)^_jO){K ^*'^\{sb)^_jB ) 


+ a^ippKuu + dd + ss)y^jO){K-^*>\{sb)y_jB-) 

+ (04 + aio){K-^*^\{su)y_jO){pp\{ub)y^jB-) 

- 2(06 + a8)(^“^*7(l + 75)M|0)(pp|h(l - -f 5 )b\B~) 
+ (04 + aio){K-^*'^pp\{su)y_jO){ 0 \{ub)y^jB~) 

- 2(06 + a 8 ){K~^*^pp\s{l + 75)M|0)(0|h(l - 75)6|S") 


ag{pp\{uu — -dd 


2 ' 5 ' 5 J v+A I 


)v*JO)(/rW|(a),.JB') 


(4.43) 


pff 

0 - 2 i = Co,; + 


Nr_ 


-c 


.eff 

2i-l) 


a2i-l — C2f_i + 




.eff 

^2i ■ 


(4.44) 


T=rO 


In Eq. ( |4.43|) the last two terms correspond to weak annihilation. As in the decay B 
pn'K~, since we do not know how to evaluate the 3-body hadronic matrix element, we will 
instead evaluate the corresponding low-lying pole diagrams with the strong process B~ 
{A[*\ followed by the weak decays K~^*^p [cf. Figs. 2(a) and 2(c)]. 


Consider the 4 intermediate poles and the hnal state K Erst. Applying factorization to 
Af, —> K~p yields 


{K-p\nw\Ab) = 


Gf 

72 ^ 


VubV*,ai - VtbV;^{a4 + aio) {K |(sM)^_^|0)(p|(M6)^_^|Afc) 


2V)ftf^*(a6 + Og) 




{K |(sm)^^^|0)(p|(m6)^^^|A6) 


rubrns 

where we have applied equations of motion 

- id^{qi^f,q2) = {mi - m2)qiq2, {^ 111 , 1 ^ 2 ) = {mi + m2)gi75g2- 


(4.45) 


(4.46) 


The pole amplitude then has the form 
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A{B —> ppK 


= 9ao^b-p Ik Uplf^^^{m]^)h[2pK ■ Pp + ^ximA, - mp) 


75 


2pK ■ Pp - ^ximA, + rUp) 

+ (A5-S^), 


>Vp X 


(Pp+Pk)^ 

(4.47) 


with 


h — VubV* tti — VtbV^*^a4 + Oio + 2(a6 + as)— 

'' rnbrus 

2 

h' = VubV* ai — VtbV^*Aa4 + aio — 2(a6 + as) — 


(4.48) 


Since pA^^s-p = g-so^B-p [Eq- (|3.23|) ], it is evident that the pole contributions arising 
from the Ab and intermediate states are constructive and dominated by the former one. 

The amplitude of B~ ppK*~ is similar to that of ppK~ except that there are no 
ae and as penguin contributions to h or h' given in Eq. (|4.48|) owing to the fact that 
(77 *“|sm| 0) = 0. For numerical calculations of decay rates we use the running quark masses 
mb{mb) = 4.4 GeV, ms{mb) = 90 MeV and the decay constant fx* = 221 MeV. Note that 
the corresponding running strange quark mass at /i = 1 GeV is 140 MeV. Applying the 
baryon form factors given by Eqs. (|4.28|) and (^4.29|) we obtain 


B{B- ppK-) = 4.0 X 10“® (2.3 x 10“^), 
B{B- ppK*-) = 2.3 X 10“® (2.1 x 10“^), 


(4.49) 


for n = 1 (n = 2), where use of the strong coupling Is'Ag^s-pl = 16 has been made. As 
stressed before, the large chiral enhancement of penguin contributions characterized by the 
x/ij^b'iT^s) term accounts for the sizable decay rate of B~ —>• ppK~. 


m 


An observation of this mode has recently been reported by Belle 
B{B- ppK-) = (4.3lJ;^ ± 0.5) x 10“®. 


(4.50) 


This is the first ever measurement of the penguin-dominated charmless baryonic B decay. 
Evidently, the model prediction is in good agreement with experiment provided that the 
baryon form factor dependence is of the monopole form (i.e. n = 1). However, in view 
of many assumptions and uncertainties involved in the calculation, the statement about the 
monopole dependence for heavy-to-light baryonic form factors should be regarded as a 
suggestion rather than a firm one. The absence of penguin contributions of qq and ag to 
K* production explains why the ppK*~ rate is smaller than ppK~, contrary to the case of 
B^ —>■ np7i~^{p'^) where the ratio of p'^ can be as large as 5. 

In Fig. 3 we show the differential decay rate dV/dt of B~ —> ppK~ where t = {pp + 
PpY = {Pb — PkY- Evidently, the spectrum peaks at f ~ 5.5 GeV^, indicating a threshold 
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FIG. 3. The differential decay rate of B —> ppK where t = {pp + Pp)'^ = {pb — PkY 


enhancement for baryon production and a fast recoil kaon accompanied by a baryon pair 
with low invariant mass. 

For other NNK^*'^ modes, it is easily seen that the pole amplitude of ^ pnK~^*^ is 
very similar to that of B~ ppK~^*^ except the A;, and poles contribute destructively 
owing to the relation q, = —Sy/Sq^n [Eq. (|3.201)1. Repeating the same calculation 

as before gives 


B{B 

B(^ 


pnK ) = 1.9 X 10 ® 
pfiK*-) = 1.8 X 10“® 


(1.5 X 10“^) 
(1.9 X 10“^) 


(4.51) 


for n = 1 (n = 2). As for —>■ nnK^^*\ its pole amplitude is the same as ^ pnK~^*^ 
except that the electroweak parameters Og and Oio in Eq. ([4.48|) are replaced by —and 
— |aio, respectively. Since these parameters are very small, the mode nnK^^*'^ has a similar 
rate as pnK~^*\ 

As for the decays ^ ^ ppK^^*^ and B~ nnK~^*\ their branching ratios are sup¬ 

pressed, of order a few times of 10“^ for K production and 5 x 10“® for K*. This is attributed 
to the fact that only the pole contributes and its coupling with the B meson and the 


nucleon is smaller compared to A;,. The current limit is B{B^ — ppK^) < 7.2 x 10 ® [^ . 


2 . 


Ap7r+(p+) 


This decay receives contributions from Figs. 2(b), 2(d), 2(e) and 2(g). The factorizable 
amplitude from Figs. (2b) and 2(d) including tree and penguin transitions is 
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(4.52) 


A{B Ap 7 r+(p+))fact = 


G^{ 


V2\ 

VthVZ 


VubV*,ai (7r+(p+) I iub)^_^\B°) {Ap\ |0) 

04 + aio)(7r+(p+)|(M6)^_^|)^)(Ap|(sM),,_^|0) 


- 2(06 + 08 ){vr+(p+)|M(l - 75)6|5 )(Ap|s(l + 75 )m|0) 


+ (04 - ^aio){7r^{p^)Ap\{sd)^_jO){0\{db)^_jB^) 

- 2(o 6 - ^08)(7r+(p+)Ap|s(l + 75)ci|0)(0|J(l - 75)&|5°) 


where the hrst term corresponds to external Ph-emission, second and third terms to the 6 —»• s 
pengnin transition and the last two terms to pengnin-indnced weak annihilation. We shall 
neglect the weak-annihilation contribntions in the practical calcnlation. Applying eqnations 
of motion we obtain 

(Ap|s(l + 75 )m| 0) = %-^-^^(Ap|s7^6|0) + (Ap|s7^75&|0) 


m.o 


m„ 


Tfl/^ Trip „Ap 


m.. 


m, 


/i ^{t)uAVp + 


rus + rrip 
1 


TUs + rrii, 


(roA + mp)pf^(t) 


+ 


t 


dm 


rriA + rup 




(4.53) 


where t = (pA+PpY and we have taken the SU(3) symmetry resnlt = 0 [see Eqs. (|4.5| ) 

and (^4.7|) ]. Since the psendoscalar form factor corresponds to a kaon pole contribntion to 
the Ap axial matrix element, it follows that 


dm = 


t — m 


(4.54) 


K 


Consequently, 


/A-i-n , ^ In\ Ap/,5^ rriA + mp ttik 

Ap S 1 + 75 M 0 = -^/i ^ t WAt^p-^^ : 

, rUs + mu t — m‘ 


pf^(^)MA75^'p• (4.55) 






K 


It is easily seen that the hrst term on the right hand side satishes the relation of vector 
current conservation in the SU(3) limit, while the second term respects the PCAC relation. 
Therefore, the above expression has a smooth chiral behavior in the zero light quark mass 
limit mg, mu 0. Applying equations of motion again yields 


(7r+|M(l - -f5)b\B ) = 


m 


B 


mz 




mb 

mr. 


(p+|m(1 - -ib)b\B ) = 2i ■ p^), 

mb 


(4.56) 


where use of Eqs. ( [f.l71 - |TT9|) has been made. Therefore, the third term in Eq. ([4.52|) is 
reduced to 
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( 7 r+|M(l - 75 ) 61 #)(Ap|s(l + 75 )m| 0 ) = 


m 


B 


m: 




mb 


Tflp^ Trip 


■- m. — m,. 


mf^ + mp mi ^ 
-W75 


ms + mu t — ml 
,m 


K 




{p^\u{l - 75 ) 6|5 )(Ap|s(l + 75 )m| 0 ) = 2i—Ao^{t){e* ■ Pb)ua 

mb 


mA TBlp yAp/^\ 

'-ms- mu ^ 


mA + mp mi Ap/,x 

- —^ T- t9i itn5 

ms + mu t - mi 

The Clebsch-Gordon coefficients for weak Ap form factors are 


Vp. 


(4.57) 


Vg 


f Ap _A 

^ v/6' 


Hence, 


(4.58) 




fAt) = -AFm 


(4.59) 


In the large t regime, the dominated axial form factor is 

/3 




01,(1). 


(4.60) 


As the B —>■ np7i^{p~^) decays, we consider two distinct empirical fits for the proton magnetic 
form factors denoted by CHT and BKLH. Using the same rnnning qnark masses as before 
we show the results of branching ratios in Table IV with and without the contributions from 
the axial form factor gi^. When including the contribution from axial form factors we shall 
assume the validity of the relation (|4.6CI|) for all the range of t. We see that the predictions 
are quite sensitive to the baryonic form factors and It is evident from Table IV 
that the factorizable contributions to # ^ Ap 7 r+(p+) are generally smaller than 1 x 10 “®, 
while the current limit is 1.3 x 10“® [^. Thus far we have neglected the nonfactorizable 
contributions from Fig. 2(e). The corresponding pole diagrams involve and S’*" poles. 
Unfortunately, it is not easy to evaluate the nonfactorizable weak matrix elements. It is 
conceivable that the total decay rate will be enhanced by a factor of 2 . At any rate, we 
conclude that the branching ratios of B ^ Apir^ are at most on the verge of 10“®. 

In Sec. IV.C below we shall explain why this penguin-dominated decay does not have 
a large rate. In contrast, the radiative baryonic decay B~ —> Apy is likely to have an 
appreciable decay rate for two reasons. First, the main pole diagram for this radiative 
decay comes from the strong process B~ AbP followed by the weak radiative transition 
Ab —>■ Ay. Since the latter is induced by the electromagnetic penguin mechanism b —>■ sy, it 
has a magnitude of order 1 x 10“® [^]. Second, the coupling of the Ab with the B meson 
and the nucleon is large. Our study indicates that B{B~ —> Apy) (1 ~ 5) x 10“® [p 6 | . 
Therefore, experimentally it would be quite interesting to measure the radiative baryonic B 
state Apy and compare with Ap7r^{p'^). 
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TABLE IV. Same as Table III except for ^ —> Ap7T~^{p'^). The axial form factor gi^{t) is 
taken to be the asymptotic form given by Eq. (4.60). 



G^m 

(CHT) 

rip 

'^M 

(BKLH) 

MS 

BSW 

MS 

BSW 

B - 

-> ApTr"*" 

2.2 X lO"’^ 

3.4 X 10“'^ 

8.0 X 10“® 

1.2 X 10“'^ 



2.9 X 10"^ 

4.3 X 10“'^ 

8.5 X 10“® 

1.3 X 10“'^ 

:b° - 

-> App+ 

2.3 X 10“^ 

3.3 X 10-'^ 

6.6 X 10“® 

9.2 X 10“® 



4.8 X 10"^ 

6.4 X 10“'^ 

8.5 X 10“® 

1.2 X 10“’' 


3. E^p7r+{p+) and 5° ^ E-n7r+(p+) 


There are several other interesting pengnin-dominated modes as listed in ( [4.3|) , for ex¬ 
ample S~A, S“S°}7r’''(p’'“). The calculations are very similar to that of 






■n,5 A, 

The relevant Clebsch-Gordon coefficients for weak form factors are 



1 

72’ 

A- = 

1 

“72’ 

" = 1, 

7'" = 

-1, 


d^-^ = 

1 

“To’ 

fE-A^ 

3 

To’ 

V2’ 

0 

w 

1 

[I] 

1 

72' 

(4.61) 


Hence, 


& = + 2 K 2 ), 

/iT = ^(-f’i;2+n"2). 


/iL" = -(K2 + 

fi.A° = 4(^2 - fu). 


(4.62) 


The dominated axial form factors in the large t regime are 

- ^(G'„ + 6G^), ^ I(G'„ + 6GI). 


(4.63) 


Obviously, A{B^ S“h7r+(p+)) = \/2A(B^ E°p7r+(p+)). 

It turns out that the branching ratios of b'^ —>■ H“A(S°)7r’'', being of order 5 x 10“®, 
are even smaller than the Ap7T~^ final state. Therefore, only the results for IlNn~^{p'^) are 
shown in Table V. We see that (i) branching ratios of E“n7r’''(p’'') he in the ranges of (1.0 ~ 
2.2) X 10“® and (0.6 ~ 1.6) x 10“®, respectively. Thus the ratio of p~^/tt~^ is not greater than 
unity, contrary to the case of ^ —>■ pnn~{p~). (ii) The decay rate of S“rl7r’''(p+) is two times 
as large as that of TF'pTT^{p'^), but the latter will be more easy to detect experimentally. 


4. ^ rj'Ap 
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TABLE V. Same as Table III except for B 
to be the asymptotic form given by Eq. (4.63|). 


(/?■*■). The axial form factor gi{t) is taken 



G^m 

(CHT) 

G^m 

(BKLH) 

MS 

BSW 

MS 

BSW 

b'" s 

1.0 X 10“® 

1.6 X 10“® 

1.4 X 10“® 

2.0 X 10“® 


1.1 X 10“® 

1.8 X 10“® 

1.2 X 10“® 

2.2 X 10“® 

^ s 

6.9 X 10“^ 

6.0 X 10“'^ 

1.0 X 10“® 

1.0 X 10“® 


1.2 X 10“® 

1.0 X 10“® 

1.6 X 10“® 

1.5 X 10“® 

— > S“n7r^ 

1.9 X 10“® 

3.3 X 10“® 

2.4 X 10“® 

4.1 X 10“® 


2.2 X 10“® 

3.7 X 10“® 

2.7 X 10“® 

4.5 X 10“® 

^ T,~np'^ 

1.4 X 10“® 

1.2 X 10“® 

2.0 X 10“® 

2.0 X 10“® 


2.4 X 10“® 

2.1 X 10“® 

3.2 X 10“® 

3.0 X 10“® 


It has been argued in |^] that B —> rj'BgB could be the most promising charmless 
baryonic modes; they may be comparable to the r]'K and a crude estimate yields r(i? —> 
rj'Ap) ^ 0.3r(S — rj'K). Of course, the study of rj'Ap is much more complicated than 
rj'K'. The factorizable amplitudes for the former involves several 3-body matrix elements 
that are difficult to evaluate. Another complication is that what is the role played by the 
gluon anomaly is still controversial and not clear even for r]'K modes, not mentioning the 
three-body one, p'BgB. A detailed study of 5 — rj'Ap will be presented elsewhere. 


C. When do we have T{B B 1 B 2 M) > T{B B 1 B 2 ) ? 

As discussed in the beginning of this section, the question of why some of three-body 
baryonic B decays in which baryon-antibaryon pair production is accompanied by a meson 
have larger rates than their two-body counterparts can be qualitatively understood in terms 
of the Dalitz plot analysis which indicates that, for example, the diquark ud has a very 
large invariant mass due to the V — A nature of the b udu process [^,^. If the ud 
forms a nucleon, then it will tend to form a highly excited baryon and will be seen as 
Nn'nin > 1). This explains why NN hnal states have small rates, why pA and SA states 
are leading tree-dominated and penguin-dominated two-body baryonic B decay modes, and 
why the three-body mode NNtt{p) is favored over the two-body one. From the calculations 
in Sections III and IV, we can give a more quantitative statement. 

The experimental fact that the penguin-dominated decay B~ ppK~ has a magnitude 
larger than the two-body counterpart B^ —> pp can be easily explained in the language of 
the pole model. The intermediate pole states are A[*^ and for the above-mentioned 
three-body hnal state and for the two-body one. First, the E;, propagator in the pole 
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amplitude for the latter is of order l/(mf — ml), while the invariant mass of the {pK~) 
system can be large enough in the former decay so that the propagator in the pole diagram 
is no longer subject to the same 1/ml suppression. Second, Ai, (and the anti-triplet bottom 
baryon Sf,) has a much larger coupling to the B meson and the light octet baryon B than 
Sfc [see Eq. (|3.23| )]. These two effects will overcome the three-body phase space suppres¬ 
sion to render the three-body mode dominant. The other examples in this category are 
r(i? —»• pn7r~) > r{B~ —>■ np) as shown before and r(i?~ —>■ Acpn~) > r{B —>■ Acp) as 
discussed in [^. We have shown before that r(B^ —>■ upn^) > r{B~ —> np) even though 
the pole diagram for the former does not have a A;, pole. This can be comprehended from 
the observation that the former is dominated by the external IT-emission contribution gov¬ 
erned by the parameter oi, while the latter proceeds via the internal W emission process. 
If the aforementioned conditions are not satisfied, then the three-body mode will not nec¬ 
essarily have larger branching ratios than the corresponding two-body ones. For example, 
the penguin-dominated decays ^ ^ ppK^, nfiK^ proceed through the pole only and 

hence their rates are suppressed. The penguin-dominated decays B Ap7r+(p+) are also 
suppressed relative to ppK^*^ modes due to the lack of A^ poles. Indeed, we found their 
magnitude does not exceed 1 x 10“®. 

V. DISCUSSIONS AND CONCLUSIONS 

We have presented a systematical study of two-body and three-body charmless baryonic 
B decays. We hrst draw some conclusions from our analysis and then proceed to discuss the 
sources of theoretical uncertainties during the course of calculation. 

1. The two-body baryonic B decay B —>■ B 1 B 2 receives main contributions from the in¬ 
ternal W-emission diagram for tree-dominated modes and the penguin diagram for 
penguin-dominated processes. We evaluate the corresponding pole diagrams to calcu¬ 
late the nonfactorizable contributions. The parity-conserving baryon matrix elements 
are estimated using the MIT bag model. We found that the bag-model estimate of 
baryon matrix elements are about three times as small as the previous calculation 
based on the harmonic oscillator model. The predicted branching ratios for two-body 
modes are in general very small, typically less than 10“®, except for the case with a A 
resonance in the hnal state. Physically, this is because the diquark system in b decay 
has a very large invariant mass and hence it tends to form a highly excited baryon 
state such as the A and will be seen as Nmi{n > 1), for example. This also explains 
the non-observation of the NN hnal states. We found that the tree-dominated decay 
B~ pA can be of order 10“® due to the large coupling of the A with the B me¬ 
son and the octet baryon. This charmless two-body baryonic mode should be readily 
accessible by B factories BaBar and Belle. 
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2. Owing to large theoretical uncertainties with parity-violating matrix elements, we focus 
only on the parity-conserving contributions for two-body hnal states. Nevertheless, 
B~ —> np, B NA and SA are purely parity-conserving, whereas B —AA is 
purely parity-violating, provided that the quark pair is created from the vacuum with 
vacuum quantum numbers {^Pq model). These features can be tested by measuring 
decay asymmetries or longitudinal polarizations. 

3. Although three-body modes in general receive factorizable contributions, not all of 
them are calculable in practice due mainly to the lack of information for three-body 
hadronic matrix elements. Therefore, in many cases we still have to reply on the pole 
approximation to evaluate the factorizable amplitudes. 

4. For three-body modes we focus on octet baryon hnal states. The tree-dominated modes 
B^ — np7i^{p'^) have a branching ratio of order (1 ~ 4) x 10“® for the tt’*' production 
and (3 ~ 5) X 10“® for the p'^ production. Moreover, B(B^ pn7r~) ~ 3 x 10“® and 
B(B^ — pnp~) ~ 8 X 10“® are predicted. There are some theoretical uncertainties for 
the prediction of B~ ppn~ and it is conjectured to have a branching ratio of order 

10-F 


5. Assuming a monopole dependence for heavy-to-light baryon form factors, we predict 
that B{B~ ppK~) ~ 4 X 10“® and the other penguin-dominated decays B~ 
ppK*~, ^ pnK~ and B^ — pnK*~ all have the branching ratio of order 2 x 10“® 
and their NN mass spectra peak at low mass. The hrst one is consistent with the 
recent measurement of B~ ppK~ by Belle. Therefore, several B — NNK^*'^ decays 
should be easily seen by B factories at the present level of sensitivity. The study of 
the differential decay rate of B~ — ppK~ clearly indicates a threshold baryon pair 
production and a fast recoil meson accompanied by a low mass baryon pair. 

6. The predictions of tree-dominated decays B pp/np, B NA and penguin- 
dominated modes B Sp, SA in the QCD sum-rule approach and the diquark model 
are quite different from the present work. Measurements of the above-mentioned modes 
can differentiate between the different approaches. 

7. The factorizable contributions to the penguin-dominated decays containing a strange 
baryon, e.g., B^ Il^p7r^{p~^), J]~n7r~^{p'^), Ap7r+(p+), are calculable. While the 
SAtt’*' state has a sizable rate, of order (1 — 3) x 10“®, the branching ratios of B ^ 
Ap7T^{p'^) are in general smaller than 10“®. 

8. Some of charmless three-body hnal states have a larger rate than their two-body coun¬ 

terparts because (i) the propagator in the pole diagrams for the three-body hnal state 
is not suppressed by f/ml, and (ii) in general the pole diagram of the former contains 
a Af, or Sfe intermediate state which has a large coupling to the B meson and the light 
baryon, for example > pn7i~) > T{B~ —> np), T{B~ ppK~) > r()^ pp), 
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or (iii) some three-body baryonic decays are dominated by the factorizable external W- 
emission governed by the parameter oi, for example, np7r~^) > r(B~ np). 

Needless to say, the calculation of baryonic B decays is rather complicated and very much 
involved and hence it suffers from several possible theoretical uncertainties. Though most of 
them have been discussed before, it is useful to make a short summary here: (i) Since it is 
very difficult to evaluate nonfactorizable and even some of factorizable amplitudes, we have 
relied on the pole approximation that, at the hadron level, these amplitudes are manifested as 
the pole diagrams with low-lying one-particle intermediate states. We use the bag model to 
evaluate the baryon matrix elements. Owing to the technical difficulties and the unreliability 
of the model for describing negative parity resonances, we limit ourselves to poles and 
hence consider only parity-conserving amplitudes. In the future we need a more sophisticated 
method to evaluate both PC and PV weak baryon matrix elements. Another important 
issue is that the intermediate pole state may be far from its mass shell and this will affect 
the applicability of the quark-model estimate of baryonic matrix elements, (ii) We have 
applied the quark-pair-creation model to estimate relative strong coupling strengths. 
This amounts to treating the strong BBi,B coupling as point-like or assuming its relative 
magnitude not being affected by the momentum dependence. However, it is not clear to us 
how good this approximation is. In the future, it is important to have a solid pQCD analysis 
to understand this issue, (iii) Heavy-to-light baryon form factors are evaluated in the non- 
relativistic quark model at zero recoil. However, their dependence is basically unknown. 
We have resorted to the pole dominance approximation by assuming a simple monopole 
or dipole momentum dependence. The unknown momentum dependence for baryon form 
factors is one of the major theoretical uncertainties, (iv) We have applied SU(3) symmetry to 
relate the octet-octet baryonic vector form factors to the magnetic and electric form factors 
of the nucleon. Experimentally, one certainly needs measurements of nucleon (especially 
neutron) electromagnetic form factors for a large range of q^. Theoretically, it is important 
to know how important the SU(3) breaking effect is and how to treat the baryonic axial 
form factors, (v) The three-body decays usually proceed through several quark diagrams. 
To simplify the calculation and to catch the main physics, we have often focused only on the 
leading factorizable quark diagrams. It remains to investigate nonfactorizable contributions 
to see their relevance. 

To conclude, we have pointed out several promising charmless two-body and three-body 
baryonic B decay modes which have branching ratios in the range of 10“® ~ 10“® and hence 
should be measurable by B factories. 
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APPENDICES 


A. BARYON WAVE FUNCTIONS 

We list the spin-flavor wave functions of baryons relevant for our purposes: 

M “ bdu)xA + (12) + (13)], 

“ bsu)xA + (12) + (13)], 

+ bsu)xs + (12) + (13)], 

= -^[buuxs + (12) + (13)], 

= -^[(bud + bdu)xs + (12) + (13)], 

Sf, ^ = -^\bddxs + (12) -I- (13)], 

S+^ = -^[suuxs + (12) + (13)], 

= -^[{sud + sdu)xs + (12) + (13)], (Al) 

S ^ = -^[sddxs + (12) -I- (13)], 

= -^[{sud - sdu)xA + (12) + (13)], 
v6 

= -^[duuxs + (12) + (13)], 

= -^[uddxs + (12) + (13)], 

A++^ = UUUX-, 

= -^[duux + (12) + (13)], 

A°^ = -^j^uddx + (12) -|- (13)], 

where abcxs = (2al6^c^ — a)Wc^ — a}h^c^)/abcxA = {cdb^'c^ — a)h^c^)/\/2^ abcx = 
(alfe^c^ + a)h^c^ + a^6lc^)/v^, and {ij) means permutation for the quark in place i with 
the quark in place j. The spin-flavor wave function of the A is expressed for = 1. The 
relative sign of baryon-pseudoscalar couplings is then fixed. 
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B. BARYON MATRIX ELEMENTS IN THE BAG MODEL 


Some of the details for evaluating the baryon matrix elements in the MIT bag model are 
already shown in Here we add the result for the matrix element of {V — A){V + A) 

current. Consider the four-quark operator O = {qa(lb)v-Ai^cQd)v-A- t)e written as 

0{x) = G{qaqb)l,_^{qcqd)y_^ where the superscript i on the r.h.s. of O indicates that the 
quark operator acts only on the ith quark in the baryon wave function. Applying the 
relations 

{qI^oIq) = UU v'v, 

{q\Ao\q) = —i{uv — v'u)a ■ f, 

{q\V\q) = —{uv -f v'u)a x f — i{uv — vu)f, (Bl) 

{q'\A\q) = {uu — v'v)a + 2v'v fa ■ f, 

leads to the PC matrix elements 

Jr^dr{q^q^\{qaqbYy_^{qcqd)l_JqU^)pc = Xi + X2 + (Xi - X2)ai ■ Us - 2 Xi(ai ■ f)(a2 ■ r) 

= Xi + X 2 + ^(Xi - 3X2)ai ■ a2, (B2) 


for {y — A){y — A) current, where we have used the relation 

5,; 


dVtfifj = — / dO, 


(B3) 


and Xi, X 2 are the bag integrals 


rR 


Xi = / r^dr[ua{r)vb{r) - Va{r)ub{r)][uc{r)vd{r) - Vc{r)ud{r)], 


rR 


X 2 = / ryr[ua{r)ub{r) + Va{r)vb{r)] [uc{r)ud{r) + Vc{r)vd{r)], 


(B4) 


with Uq{r) and Vq{r) being the large and small components, respectively, of the lS' 1/2 quark 
wave function (see |p8|,|^ for detail). Likewise, for {V — A){V + A) current we obtain 

Jr^dr{qlqyqaqb)l_^{qcqd)lyjq'lq 2 )pc = {Xi + X2) + ^[Xi + X2 - 2 (X( - X^di ■ 0^2, (B 5 ) 


with 


rR 


X[= r^dr[ua{r)vb{r) +Va{r)ub{r)][uc{r)vd{r) +Vc{r)ud{r)], 


rR 

X' = / r^dr[ua{r)ub{r) - Va{r)vb{r)][uc{r)ud{r) - Vc{r)vd{r)]. 

Jo 

For numerical estimates of the bag integrals, we shall use the bag parameters 

RT'u = RT-d = 0, ms = 0.279 GeV, me = 1.551 GeV, mb = 5.0 GeV, 

= 2.043. X, = 2.488. = 2.948. 14 = 3.079, R = 5.0 GeV^'. 


(B6) 


(B7) 
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